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ABSTRACT 

We  use  high  dynamic  range,  high-resolution  L-band  spectroscopy  to  measure  the  radial  velocity  (RV)  variations  of 
the  hot  Jupiter  in  the  r  Bootis  planetary  system.  The  detection  of  an  exoplanet  by  the  shift  in  the  stellar  spectrum 
alone  provides  a  measure  of  the  planet’s  minimum  mass,  with  the  true  mass  degenerate  with  the  unknown  orbital 
inclination.  Treating  the  r  Boo  system  as  a  high  flux  ratio  double -lined  spectroscopic  binary  permits  the  direct 
measurement  of  the  planet’s  true  mass  as  well  as  its  atmospheric  properties.  After  removing  telluric  absorption  and 
cross-correlating  with  a  model  planetary  spectrum  dominated  by  water  opacity,  we  measure  a  6a  detection  of  the 
planet  at  Kp  =  111  ±  5  km  s'1 ,  with  a  la  upper  limit  on  the  spectroscopic  flux  ratio  of  1CT4.  This  RV  leads  to 
a  planetary  orbital  inclination  of  i  =  45434°  and  a  mass  of  Mp  =  5.904q350  Mjup.  We  report  the  first  detection  of 
water  vapor  in  the  atmosphere  of  a  non-transiting  hot  Jupiter,  r  Boo  b. 
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1.  INTRODUCTION 

Since  the  first  detection  of  an  exoplanet  around  a  main- 
sequence  star  (Mayor  &  Queloz  1995),  astronomers  have  dis¬ 
covered  hundreds  of  exoplanets  using  the  radial  velocity  (RV) 
technique  (Wright  et  al.  2012).  This  powerful  tool  for  discov¬ 
ering  exoplanets  only  measures  the  minimum,  or  “indicative” 
mass,  MpSini,  leaving  a  degeneracy  between  two  interesting 
properties  of  the  system,  one  physical  and  one  orbital.  While  the 
indicative  mass  is  useful  for  statistical  studies,  uniquely  measur¬ 
ing  the  true  mass  of  an  exoplanet  not  only  yields  a  key  physical 
property,  but  also  furthers  our  understanding  of  planetary  forma¬ 
tion  and  evolution  via  measurement  of  the  true  mass  distribution 
of  exoplanets  (Zucker  &  Mazeh  2001 ;  Weiss  et  al.  2013).  For  ex¬ 
ample,  Ho  &  Turner  (201 1)  demonstrate  that  knowledge  of  the 
true  mass  distribution  is  necessary  to  convert  a  minimum  mass, 
Mp  sin  i,  into  an  estimate  of  a  planet’s  true  mass  for  RV-detected 
systems. 

One  class  of  objects  for  which  the  planet  mass  can  be  directly 
determined  are  those  that  transit  their  host  star.  Hundreds 
of  transiting  planets  have  been  discovered  and  characterized, 
and  the  ongoing  Kepler  mission  has  found  potential  exoplanet 
candidates  numbering  in  the  thousands  (Borucki  et  al.  2011; 
Batalha  et  al.  2013).  Transit  events  also  provide  atmospheric 
information  through  transmission  spectroscopy  and  secondary 
eclipses.  Investigators  have  measured  spectra  of  some  of  the 
larger  transiting  planets  (Charbonneau  et  al.  2002),  leading  to 
the  discovery  of  species  such  as  water,  methane,  and  carbon 
monoxide  (Knutson  et  al.  2012;  Berta  et  al.  2012;  Crossfield 
et  al.  2013;  Baskin  et  al.  2013).  A  variety  of  spectral  retrieval 
methods  have  been  used  to  verify  the  results,  confirming  the 
hypothesis  that  O-  and  C-bearing  gases  are  present  in  the 
atmospheres  of  hot  Jupiters  (Madhusudhan  &  Seager  2009;  Line 
et  al.  2012). 

Recently,  a  technique  previously  used  to  detect  low  mass 
ratio  spectroscopic  binary  stars  has  been  applied  to  stars  known 


to  host  exoplanet  systems.  The  direct  RV  detection  of  an  exo¬ 
planet  involves  separating  the  planetary  and  stellar  components 
spectroscopically.  The  high  flux  ratio  between  the  primary  and 
companion  makes  these  detections  difficult,  but  not  impossi¬ 
ble,  thanks  to  modern  infrared  echelle  spectrographs.  Snellen 
et  al.  (2010)  detected  CO  on  HD  209458  b  with  a  precision  of 
2  km  s_1  using  the  R  =  100,000  Very  Large  Telescope  CRIRES 
instrument.  This  detection  provided  the  spectroscopic  orbit  of 
the  system  and  determined  the  true  mass  of  the  planet. 

Indeed,  this  technique  can  be  applied  to  non-transiting, 
RV-detected  exoplanets  to  extract  the  unknown  inclination  and 
true  mass.  CRIRES  was  also  used  to  detect  CO  on  r  Bootis  b 
(Brogi  et  al.  2012),  the  first  ground-based  detection  of  a  short- 
period  non-transiting  exoplanet  atmosphere,  a  result  confirmed 
shortly  thereafter  by  Rodler  et  al.  (2012).  These  studies  provide 
the  true  mass  of  the  planet  and  probe  the  chemical  composition 
of  its  atmosphere.  A  combination  of  high  signal-to-noise,  high 
spectral  resolution,  and  coverage  of  multiple  CO  overtone 
lines  was  required  to  achieve  the  sensitivity  required  for  these 
detections. 

Despite  the  agreement  between  the  two  groups,  the  direct 
detection  of  exoplanets,  especially  r  Boo  b,  has  a  storied 
history.  Collier  Cameron  et  al.  (1999)  first  reported  a  detection 
of  reflected  light  from  r  Boo  b  more  than  a  decade  ago  and 
Wiedemann  et  al.  (2001)  reported  a  possible  detection  of  CH4 
in  the  planet’s  atmosphere  soon  thereafter,  using  the  same  orbital 
solution.  However,  the  planetary  velocity  from  the  earlier  results 
disagrees  with  more  recent  findings  discussed  above. 

Here  we  report  a  detection  of  water  vapor  in  the  atmosphere 
of  r  Boo  b,  using  spectroscopic  observations  centered  around 
3.3  (im.  The  r  Bootis  system  is  a  stellar  binary  comprised  of 
an  F-type  “A”  component  and  an  M-dwarf  “B”  component. 
The  hot  Jupiter  r  Boo  b  orbits  the  larger  “A”  component 
with  a  period  of  3.312  days  and  a  minimum  (indicative)  mass 
M  sin  i  =  3.87  Mjup  (Butler  et  al.  1997).  Our  detection  confirms 
the  mass  determined  by  Brogi  et  al.  (2012)  and  serves  to 
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Table  1 

/,-band  Observations  of  Tau  Boo  b 


Date 

JD  2450000 

Phase 

(rad) 

^bary 

(km  s-1) 

S/N3^m 

2011  May  21 

5702.8542 

0.3199 

-17.36 

6525 

2012  Apr  3 

6021.0625 

0.3847 

2.26 

4768 

2012  Apr  1 

6019.0625 

0.7810 

3.02 

2871 

2011  Mar  14 

5634.9375 

0.8164 

11.75 

5467 

2011  Mar  24 

5644.9375 

0.8353 

7.39 

6575 

further  characterize  the  atmospheric  chemistry  of  this  exoplanet. 
Data  from  five  epochs  reveal  orbital  motion  of  r  Boo  b  that 
is  consistent  with  that  reported  by  Brogi  et  al.  (2012).  The 
planetary  template  spectrum  used  in  the  cross-correlation  is 
dominated  by  water  vapor  opacity,  providing  strong  evidence  of 
water  in  the  atmosphere  of  a  non-transiting  hot  Jupiter  for  the 
first  time. 

2.  METHODS 

2.1.  Observations  and  Data  Reductions 

We  observed  the  r  Boo  system  on  five  separate  nights  in 
2011  March,  2011  May,  and  2012  April  (Table  1),  chosen 
to  optimize  phasing  near  orbital  quadrature  and  allowing  for 
maximum  separation  of  stellar  and  planetary  lines.  We  use 
the  Near  Infrared  Echelle  Spectrograph  (NIRSPEC;  McLean 
et  al.  1995)  at  the  W.  M.  Keck  Observatory,  which  provides 
high  resolution  ( R  ~  25,000  for  a  three-pixel  slit)  in  multiple 
orders  at  the  wavelengths  of  interest  to  study  near-infrared 
water  emission.  We  obtained  spectra  covering  3.404-3.457  //m, 
3.256-3.307  jim,  3.121-3.170 /xm,  and  2.997-3.043  /xm.  Each 
epoch  covers  a  total  elapsed  time  of  approximately  1  hr  and  is 
comprised  of  a  continuous  sequence  of  hundreds  of  exposures. 

We  extract  NIRSPEC  spectra  using  a  custom  Interactive  Data 
Language  optimal  extraction  pipeline,  similar  to  that  described 
by  Cushing  et  al.  (2004).  The  spatial  profile  weighting  intrinsic 
to  optimal  extraction  provides  a  reliable  method  for  detecting 
and  removing  bad  pixels  due  to  detector  defects  and  cosmic 
rays.  It  adjusts  for  seeing  variations  that  occur  over  the  course 
of  an  observation,  and  can  also  minimize  the  contamination  from 
nearby  stars  that  happen  to  fall  on  the  slit.  We  have  determined 
that  the  contamination  from  the  M-dwarf  companion  r  Boo  B 
is  negligible  in  all  epochs  of  our  extracted  r  Boo  A  spectra. 
Arc  lamps  typically  used  for  wavelength  calibration  do  not 
provide  useful  reference  lines  in  the  thermal  infrared.  Instead, 
we  wavelength  calibrate  our  spectra  using  unblended  telluric 
features  with  accurately  known  rest  wavelengths.  Between  15 
and  30  individual  telluric  lines  are  used  for  each  order. 

We  first  correct  the  bulk  telluric  absorption  with 
TERRASPEC  (Bender  et  al.  2012),  a  synthetic  forward¬ 
modeling  algorithm  that  uses  the  line-by-line  radiative  transfer 
model  LBLRTM  (Clough  et  al.  2005)  to  generate  a  synthetic 
absorption  function  for  the  Earth’s  atmosphere  (hereafter  TAF). 
The  TAF  includes  continuum  absorption  and  line  absorption 
for  28  molecular  species,  with  line  information  provided  by  the 
HITRAN  2008  database  (Rothman  et  al.  2009).  Vertical  mix¬ 
ing  profiles  for  the  seven  most  prominent  species  (H2O,  CO2, 
O3,  N2O,  CO,  CH4,  and  O2)  can  be  scaled  or  adjusted  from  the 
LBLRTM  default  profiles,  which  include  the  US  Standard  1976, 
tropical,  midlatitude,  and  subartic  models.  These  models  rep¬ 
resent  the  average  atmosphere  for  their  respective  latitudes  (in 
1976),  which  is  a  good  initial  guess  for  the  instantaneous  atmo¬ 
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sphere  corresponding  to  a  single  observation.  We  use  the  trop¬ 
ical  model  to  provide  initial  vertical  profiles  for  spectroscopy 
obtained  from  Mauna  Kea. 

TERRASPEC  convolves  the  TAF  with  an  instrumental  broad¬ 
ening  profile  (hereafter  IP),  which  is  measured  from  the  data. 
The  IP  is  parameterized  as  a  central  Gaussian  surrounded  by 
satellite  Gaussians  offset  by  a  fixed  percentage  of  the  Gaussian 
width,  and  with  adjustable  amplitudes.  This  is  similar  to  the  IP 
parameterization  described  by  Valenti  et  al.  (1995)  for  use  with 
h  absorption  cells.  Typically  two  to  four  satellite  Gaussians 
serve  to  model  the  IP.  The  instrumental  broadened  TAF  is  then 
multiplied  by  a  low-order  wavelength  dependent  polynomial  to 
correct  for  the  combined  effects  of  blaze  sensitivity  and  stellar 
continuum.  TERRASPEC  uses  the  least-squares  fitting  algo¬ 
rithm  MPFIT  (Markwardt  2009)  to  optimize  the  parameters  com¬ 
prising  the  TAF,  IP,  and  continuum.  The  FWHM  of  the  profile  is 
consistently  1.3  x  10~4  /xm,  which  yields  R  —  X/ AX  —  24,000. 
Spectral  regions  containing  stellar  absorption  are  excluded  from 
the  optimization  by  adaptive  masks. 

The  initial  telluric  correction  can  be  significantly  affected  by 
the  instrumental  fringing.  We  therefore  mask  strong  stellar  lines 
from  the  telluric  corrected  spectrum,  and  analyze  the  remaining 
spectral  range  with  a  Lomb-Scargle  periodogram  (Scargle 
1982)  to  measure  the  frequency  and  power  of  the  individual 
fringes  present  in  our  spectra.  Two  prominent  fringes  are  seen 
at  ~  1 .75  cirT  'and  2.18  cm-1.  We  then  calculate  the  composite 
fringe  function  as  the  product  of  the  individual  fringes,  and 
divide  it  into  the  non-telluric-corrected  spectrum.  The  defringed 
spectrum  is  then  reprocessed  with  TERRASPEC,  yielding 
the  telluric -free  and  continuum-normalized  stellar  spectrum, 
a  measurement  of  the  TAF,  and  a  parameterization  of  the 
IP.  This  process  was  often  iterated  two  to  three  times  to 
ensure  full  removal  of  fringing  and  telluric  absorption.  Regions 
with  atmospheric  transmission  ^40%  of  the  continuum  value, 
determined  from  the  TAF,  are  masked  and  excluded  from  the 
final  planet  search.  Figure  1  demonstrates  the  telluric  removal 
process. 

2.2.  Two-dimensional  Cross  Correlation 

Next  we  use  a  two-dimensional  cross  correlation  analysis, 
TODCOR  (Zucker  &  Mazeh  1994),  to  simultaneously  extract 
the  planetary  and  stellar  velocity  shifts.  We  generated  a  synthetic 
stellar  model  of  the  r  Boo  A  spectrum  using  a  recent  version 
of  the  LTE  line  analysis  code  MOOG  (Sneden  1973)  and  the 
MARCS  grid  of  stellar  atmospheres  (Gustafsson  et  al.  2008). 
The  input  linelist  was  created  by  detailed  matching  of  a  synthetic 
solar  spectrum  to  the  ATMOS  ATLAS-3  infrared  solar  spectrum 
(Abrams  et  al.  1996),  starting  from  the  solar  linelist  generated  by 
Sauval  (see  Hase  et  al.  2006).  Using  the  MARCS  solar  model 
and  solar  abundances  in  Grevesse  et  al.  (2007),  adjustments 
were  then  made  to  the  atomic  line  parameters,  in  particular  the 
^/-values  and  damping  constants,  to  fit  the  solar  spectrum. 
For  r  Boo,  a  stellar  atmosphere  with  effective  temperature 
Tefi  =  6375  K,  surface  gravity  log  g  —  4.0,  and  metallicity 
[m/H]  =  +0.25  was  adopted,  based  on  a  review  of  abundance 
analyses  in  the  literature.  Individual  abundances  were  set  by 
matching  observed  lines  for  elements  that  were  well  measured 
by  NIRSPEC  (Fe,  Si,  Mg,  Na);  otherwise,  an  abundance  of 
+0.25  was  used. 

A  plane -parallel  model  was  calculated  for  r  Boo  b  using 
the  PHOENIX  stellar  and  planetary  atmosphere  code  (Barman 
et  al.  2001,  2005).  The  planet  is  very  hot,  with  an  equilibrium 
day  side  temperature  between  1600  and  2000  K,  depending  on 
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the  observed  spectra  as  .S',  and  the  template  spectra  as  f\  = 
f(Xj  +  AA)  yields 


2  v-a  (Si  -  fi)1 


=  ±,J2(s?  +  fr-2sifi). 


(2) 


Figure  1.  Longest  wavelength  order  of  data  from  2011  March  14.  Top:  the 
original  wavelength-calibrated  data.  Middle:  telluric -removed  spectra  with 
stellar  features  overlaid  (dotted  line).  Bottom:  stellar-  and  telluric-removed 
spectrum. 


the  day-to-night  redistribution  of  incident  stellar  flux.  Only 
2i x  redistribution  (Teq  ~  2000  K)  was  used  here.  With  an 
unknown  planet  radius,  the  surface  gravity  was  arbitrarily  set  to 
104  cm  s-2.  As  discussed  above,  water  lines  are  the  primary 
signal  we  seek  in  the  planetary  spectrum  and  to  accurately 
model  these  we  use  the  best  available  water  line  list  from 
the  ExoMol  group  (Barber  et  al.  2006).  A  final  high-resolution 
spectral  template  was  calculated  at  lOx  the  observed  resolution 
(A  A.  =  0.05  A). 

For  every  epoch,  the  target  spectrum  is  cross-correlated  to 
determine  the  cross-correlation  function  (CCF)  for  each  order. 
The  planet/star  spectroscopic  flux  ratio  is  set  to  1 0  \  the 
same  order  of  magnitude  as  the  expected  photometric  contrast 
between  the  two  objects.  We  also  tested  flux  ratios  from  10-3  to 
10~7  and  the  shape  of  the  resulting  maximum  likelihood  (ML) 
function  remains  the  same  because  the  analysis  is  only  weakly 
sensitive  to  the  absolute  contrast  ratio. 


2.3.  Maximum  Likelihood  Analysis 

To  find  the  most  likely  solution,  each  CCF  must  be  converted 
into  a  probability,  or  likelihood,  £.  To  do  this,  we  start  with  a 
relationship  between  C  and  the  familiar  y2  statistic 


£  = 


\f2nai 


rM-f 


l«g£  =  *  -  y, 


(1) 


Each  part  of  the  quotient  can  be  summed  individually.  For  a 
continuum-normalized  spectrum  with  N  points,  this  results  in 

°2  =  yY.sl  <3> 

9  Sjfi 

Xi  =  const-  —■  (4) 

The  second  term  on  the  right-hand-side  of  Equation  (4)  is  simply 
the  CCF.  Thus,  the  CCF  and  £  are  related  by 

log  £  =  const  +  CCF.  (5) 

The  goal  is  to  maximize  log  £  for  both  the  stellar  and 
planetary  velocity  shifts,  where  the  sum  is  over  all  spectral 
orders  for  an  individual  epoch.  Figure  2  (top  panel)  demonstrates 
that  the  stellar  velocity  clearly  stands  out  as  the  most  likely 
solution  for  one  epoch,  as  is  the  case  for  all  other  epochs  as 
well.  Our  NIRPSEC  spectra  have  insufficient  RV  precision  to 
be  sensitive  to  the  orbital  motion  of  r  Boo  A  so  the  stellar  RV 
is  therefore  fixed  at  the  systemic  value.  Barycentric  movement 
is  accounted  for. 

The  ML  solution  of  the  planet’s  orbit  is  much  more  complex. 
The  likelihood  is  proportional  to  the  CCF,  a  two-dimensional 
surface  that  reflects  coherence  in  features  between  template 
spectra  and  those  of  the  target.  Due  to  the  multiplicity  of  ro vibra¬ 
tional  transitions  in  the  asymmetric  top  spectrum  of  hot  water, 
the  correlation  coefficient  remains  large  (>0.9)  over  significant 
offsets.  This  results  in  multiple  “peaks”  at  incorrect  velocity 
lags,  as  can  be  seen  in  Figure  2.  Thus,  a  single  epoch  would 
lead  to  a  degeneracy  of  solutions.  Fortunately,  we  do  not  worry 
about  alignment  between  planet  and  telluric  spectra  because  not 
only  are  the  spectra  very  different,  but  the  combined  systemic 
and  planetary  RV  shifts  the  planetary  spectrum  sufficiently  to 
avoid  collisions  between  spectral  features. 

2.4.  Orbital  Solution 

To  break  the  degeneracy  several  epochs  of  data  must  be 
brought  to  bear.  Using  the  known  period  of  the  planet,  we 
calculate  the  orbital  phase  of  each  epoch  (Table  1)  and  use 
that  to  find  the  most  likely  planetary  velocity  consistent  with  a 
circular  Keplerian  orbit  (the  estimated  planetary  eccentricity  is 
small,  e  ~  0.02;  Butler  et  al.  1997;  Brogi  et  al.  2012).  Since  the 
absolute  orbital  velocity  of  the  planet  is  known  from  the  period, 
we  are  actually  interested  in  Kp.  or  semi-amplitude.  A  range  of 
Kp  is  tested,  each  corresponding  to  a  different  inclination  of  the 
planet,  as  well  as  a  unique  mass.  Furthermore,  each  Kp  leads  to 
a  different  velocity  lag  at  a  given  phase,  such  that 

Up]  =  Kp  sin(&>r  +  </>)  +  y,  (6) 


where  at  is  a  constant  that  does  not  matter  when  comparing 
relative  likelihoods,  assuming  cr,-  =  o  =  const.  Denoting 
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where  a>  =  2n / P,  (p  is  a  phase  lag  that  can  be  set  to  zero  by 
choosing  the  proper  starting  date,  and  y  is  the  combined  stellar 
barycentric  and  systemic  velocities  ( V  sin  i  ~  15  km  s  '  1  (Butler 
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Velocity  Shift  (km/s) 

Figure  2.  (A)  ML  function  of  the  stellar  velocity  shift  using  the  data  from  201 1 
March  24.  (B)-(F)  ML  function  for  the  planet  signal  for  201 1  March  24,  201 1 
March  14,  2012  April  1,  201 1  May  21,  and  201 1  April  3,  respectively.  Note  the 
changing  sign  of  planetary  velocities  for  each  epoch,  which  are  anti-correlated 
with  the  sign  of  the  stellar  RV  shift.  The  vertical  lines  correspond  to  the  velocity 
shift  at  a  given  epoch  for  an  orbital  solution  with  Kp  =  70  (dashed),  Kp  =  90 
(dotted),  and  Kp  =  1 10  km  s-1  (solid). 

et  at.  2006)  and  Vbary  listed  in  Table  1).  We  then  seek  the  best- 
fitting  Kp  by  maximizing  the  sum  of  the  likelihood  of  vp\  for  all 
epochs,  whose  highest  value  indicates  the  most  likely  solution 
for  Kp. 

3.  RESULTS 

It  is  clear  that  the  ML  curve  for  each  epoch  in  Ligure  2  shows 
multiple  peaks.  Only  one  of  these  peaks  per  epoch  can  be  the 
true  peak  that  results  from  the  planet  signal  in  our  spectra; 
the  other  peaks  are  artifacts  caused  by  chance  misalignments. 
The  molecular  lines  in  the  model  water  spectrum  can  align 
by  chance  with  both  signal  and  noise  features  in  the  observed 
spectra  and  with  other  molecules  in  the  planetary  atmosphere 
that  have  not  been  included  in  the  planetary  template.  However, 
the  true  peaks  can  be  distinguished  from  artifacts  by  requiring  a 
multi-epoch  orbital  solution  that  is  consistent  with  a  Keplerian 
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Planetary  Radial  Velocity  (km/s) 

Figure  3.  Normalized  log  likelihood  as  a  function  of  planetary  velocity,  Kp. 
(A)  Results  from  synthetic  spectra,  composed  of  the  stellar  and  planetary 
templates,  with  a  planetary  signal  injected  at  70  km  s-1  and  analyzed  with 
the  same  procedures  applied  to  the  data.  (B )  Data  analyzed  using  a  planet-to- 
star  flux  ratio  of  10~5  for  a  water  vapor  around  a  planet  with  Teq  ~  2000  K. 
(C)  Same  as  (A)  but  for  a  signal  injected  at  1 10  km  s-1 . 

orbit.  The  dashed,  dotted,  and  solid  lines  in  Ligure  2  represent 
the  tipi  for  three  Keplerian  solutions  with  Kp  =  70  km  s_1, 
90  km  s_  1 ,  and  110  km  s_  1 ,  respectively.  While  none  correspond 
to  the  highest  peak  for  all  epochs,  certain  orbital  solutions  find 
peaks  more  often  than  troughs. 

To  find  the  most  likely  orbit,  the  likelihoods  at  every  epoch  for 
a  given  Kp  are  combined.  The  second  panel  in  Ligure  3  presents 
this  composite  ML,  along  with  the  corresponding  planetary 
mass.  This  function  represents  the  sum  of  the  log  likelihoods  of 
all  five  epochs.  There  are  two  peaks  in  the  likelihood  function, 
at  Kp  ~  70  km  s_1  and  111  km  s_1.  We  have  demonstrated 
empirically  the  degenerate  effects  of  cross-correlating  a  water 
spectrum.  Now  we  explore  the  systematics  of  the  analysis  and 
consider  how  both  the  properties  of  the  individual  stellar  and 
planetary  template  autocorrelation  functions,  as  well  as  the 
correlation  of  the  two  template  spectra,  can  affect  the  composite 
likelihood  function. 

To  do  this,  synthetic  data  sets  are  subjected  to  the  same 
analysis  as  the  original  data.  Synthetic  target  spectra  of  just 
the  stellar  and  planetary  templates  are  used,  each  shifted  to 
the  correct  velocity  for  a  given  epoch.  The  results  are  given  in 
Ligure  3,  whose  top  panel  shows  the  normalized  log  likelihood 
from  a  planetary  spectrum  injected  at  70  km  s'1 ,  with  a  planet- 
to-star  flux  ratio  of  1  (T 5.  The  proper  signal  is  clearly  retrieved. 
However,  the  third  panel  shows  that  at  an  injected  planetary  RV 
of  1 10  km  s-1  and  the  same  contrast  ratio,  a  double-peaked  log 
likelihood  results.  Thus,  following  exactly  the  same  procedure 
as  was  used  for  the  observations,  a  “perfect”  target  spectrum 
with  no  noise  and  no  (terrestrial)  atmosphere  retrieves  both  the 
correct  and  an  additional  signal,  with  a  structure  that  mimics 
the  data,  likely  as  a  result  of  the  complex  hot  water  spectrum 
near  3  /jm. 
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Planetary  Velocity  (km/s) 


Figure  4.  Normalized  log  likelihood  as  a  function  of  planetary  velocity,  Kp,  for  different  spectroscopic  flux  ratios.  Left:  results  from  synthetic  spectra,  composed  of 
the  stellar  and  planetary  templates,  with  a  planetary  signal  injected  at  70  km  s-1  and  analyzed  with  the  same  procedures  applied  to  the  data.  Right:  same  for  a  signal 
injected  at  110  km  s_1 .  The  signal  at  1 10  km  s-1  requires  a  much  stronger  planetary  signal  to  uniquely  determine  the  correct  orbital  velocity. 


At  a  sufficiently  high  planet-to-star  flux  ratios,  the  true  plan¬ 
etary  signal  should  dominate  the  posterior  likelihood.  Indeed, 
using  a  flux  ratio  of  10~3,  Birkby  et  al.  (2013)  detect  water 
absorption  from  another  (transiting)  hot  Jupiter,  HD  189733  b. 
The  right  panel  of  Figure  4  shows  that  when  the  flux  ratio  of  the 
planetary  signal  is  increased  to  this  level,  the  singular  correct 
velocity  is  retrieved.  We  can  constrain  the  spectroscopic  con¬ 
trast  of  r  Boo  b  relative  to  its  host  star  by  comparing  the  data 
to  these  synthetic  fits.  The  simulations  show  that  for  a  10  4, 
the  correct  110  km  s_I  solution  should  present  a  larger  ML. 
Since  our  data  do  not  demonstrate  this,  we  can  put  a  ler  upper 
limit  on  the  high  resolution  spectroscopic  flux  ratio  of  1 0  4  at 
~3.3/xm.  The  left  panel  of  Figure  4  shows  that  a  planetary  signal 
of  70  km  s'1  would  be  uniquely  retrieved  for  all  realistic  flux 
ratios. 

4.  CONCLUSIONS  AND  FUTURE  WORK 

The  detection  of  r  Boo  b  has  been  a  difficult  quest.  Contradic¬ 
tory  results  have  been  published  over  the  past  15  years  and  the 
difficulty  in  both  thoroughly  removing  the  telluric  absorption 
and  also  extracting  the  diminutive  planetary  signal  is  not  to  be 
underestimated.  Multiple  observations  at  different  wavebands 
improve  the  validity  of  the  detection,  and  facilitate  a  variety  of 
statistical  tests  to  ensure  an  accurate  measurement. 

Here,  the  validity  of  the  most  likely  solutions  for  the  RV 
of  r  Boo  b  are  explored.  Of  the  two  prominent  velocities  that 
fit  the  data,  one  is  shown  to  result  from  the  systematics  of 
the  cross-correlation  analysis  performed  on  this  high  flux  ratio 
spectroscopic  binary  system.  Interestingly,  previous  velocity 
studies  of  the  planet  in  reflected  light  (Collier  Cameron  et  al. 
1999)  and  methane  (Wiedemann  et  al.  2001)  retrieved  values 
close  to  this  artifact.  We  suggest  that  residual  water  vapor  in 
the  atmosphere,  after  the  telluric  removal,  might  have  been 
responsible  for  the  false  value  previously  reported  from  infrared 
observations.  As  methods  for  telluric  corrections  have  improved 
over  the  past  decade,  this  problem  can  be  overcome  and  the 
correct  RV  retrieved. 

Our  analysis  gives  a  6cr  detection  of  the  planet  at  Kp  —  1 1 1  ± 
5  km  s"1  for  r  Boo  b,  with  a  ler  upper-limit  on  the  3.3  /tm 
planet/star  spectroscopic  flux  ratio  of  1 0  4.  To  determine  the 
significance  of  our  detection,  we  injected  synthetic  signals  at  a 


variety  of  planetary  velocities  and  planet-to-star  spectroscopic 
flux  ratios  and  constructed  the  chi-square  surface  of  the  ML 
fits.  The  orbital  velocity  is  in  good  agreement  with  previous  RV 
amplitude  detections  via  CO  by  Brogi  et  al.  (2012)  and  Rodler 
et  al.  (2012),  and  our  analysis  reveals  the  presence  of  water  vapor 
in  the  planet’s  atmosphere.  Furthermore,  using  a  stellar  mass  of 
1.34Uqq5349  Mg  (Takeda  et  al.  2007)  and  a  stellar  velocity  semi- 
amplitude  of  0.4664  ±  0.0033  km  s-1  (Brogi  et  al.  2012),  we 
derive  a  planetary  mass  of  MP  =  5.90jlo  20  A7jup  with  an  orbital 
inclination  of  i  —  45jl340. 

The  technique  presented  here  is  in  its  nascent  stages,  and  the 
work  is  by  no  means  complete.  Additional  quantitative  char¬ 
acterization  of  the  physical  properties  (temperature,  opacity) 
and  composition  (and  thus  estimates  of  vertical  mixing)  of  r 
Boo  b’s  atmosphere  will  require  both  significant  simulations  of 
the  atmospheric  radiative  transfer  and  data  analysis  along  with 
additional  data  at  longer  and  shorter  wavelengths.  Such  work 
is  beyond  the  scope  of  this  Letter,  but  expanded  studies  of  r 
Boo  using  the  same  methods  described  above,  but  also  includ¬ 
ing  molecules  such  as  CH4  and  CO  at  longer  wavelengths  than 
examined  here,  are  underway.  Indeed,  although  the  mole  frac¬ 
tions  of  methane  should  be  insignificant  for  r  Boo  b,  we  have 
searched  for  the  molecule  using  the  methods  described  herein 
and  have  no  detection  to  report  in  our  data.  With  careful  anal¬ 
ysis  in  the  future,  the  relative  abundances  of  these  molecules 
in  the  atmospheres  of  non-transiting  exoplanets  can  be  ascer¬ 
tained.  Further  applications  of  this  technique,  using  water  and 
methane  as  template  molecules,  to  a  variety  of  additional  hot 
Jupiter  exoplanets  will  be  reported  elsewhere. 
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